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We study the excited random walk, in which a walk that is at a site that contains cookies eats 
one cookie and then hops to the right with probability p and to the left with probability q = 1 — p. 
If the walk hops onto an empty site, there is no bias. For the 1-excited walk on the half-line (one 
cookie initially at each site), the probability of first returning to the starting point at time t scales 
as 

j-(2-p) Although the average return time to the origin is infinite for all p, the walk eats, on 
average, only a finite number of cookies until this first return when p < 1/2. For the infinite line, the 
probability distribution for the 1-excited walk has an unusual anomaly at the origin. The positions of 
the leftmost and rightmost uneaten cookies can be accurately estimated by probabilistic arguments 
and their corresponding distributions have power-law singularities near the origin. The 2-excited 
walk on the infinite line exhibits peculiar features in the regime p > 3/4, where the walk is transient, 
including a mean displacement that grows as t" , with ^ > \ dependent on p, and a breakdown of 
scaling for the probability distribution of the walk. 

PACS numbers: 02.50.-r, 05.40.-a 



I. INTRODUCTION 

In this work, we study the /c-excited random walk in 
one dimension. This recently-introduced model 0, 
gives rise to a stochastic process with many unusual prop- 
erties H, IM Q ' The definition of the excited random 
walk is as follows: Initially each site of a one-dimensional 
chain contains a fixed number of cookies k (with fc > 1). 
A random walk makes nearest-neighbor hops between 
sites on the chain, with hopping probabilities that de- 
pend on whether or not the current site contains a cookie. 
When the walker hops onto a site that contains one or 
more cookies, the walker eats a cookie and, at the next 
time step, the walk hops to the right with probability 
p > 1/2 and to the left with probability q — \~p. If the 
walker hops onto an empty site, then the hopping prob- 
abilities to the left and right at the next step are both 
1/2. This bias can be viewed as the foraging strategy 
of a primitive organism. When it encounters food, the 
organism consumes the food and continues to move in 
a preferred direction. Conversely, the organism wanders 
aimlessly when it does not encounter any food. 

A basic feature of the excited random walk is the tran- 
sitory nature of the bias; that is, when all the cookies 
at a given site are eaten the subsequent motion is un- 
biased. This non-Markovian feature is responsible for 
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FIG. 1: Space-time plot of the first few steps of a 1-excited 
random walk in one dimension with reflection at the origin 
(left), and on the infinite line (right). Cookies are shown as 
large dots. Steps occurring with probability p and q are shown 
as thick and dashed lines, while steps occurring with proba- 
bility 1/2 are shown as thin lines. For the half-line system, 
the last cookie at n = 4 is first hit after 10 time steps. 



the unusual properties of the excited random walk. Re- 
lated non-Markovian constraints arise in many important 
stochastic phenomena. Perhaps the best known exam- 
ple is the self-avoiding random walk, where returns to 
previously- visited sites are strictly forbidden 0. This 
leads to an expansion, relative to the pure random walk, 
in which the root-mean-square displacement grows faster 
than t^/^ when the spatial dimension is less than 4. A 
related example is the "true" self-avoiding walk, where 
returns to previously-visited sites are inhibited by an en- 
ergetic constraint [g . The opposite type of perturbation 
occurs in the reinforced random walk '9] , where the walk 
has a preference to hop onto previously- visited sites. De- 
pending of the functional form of this reinforcement, the 
walk may, quite surprisingly, "condense" onto a support 
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that consists of only a finite number of sites as t — s- oo 
i 

The excited random walk represents a simple and sur- 
prisingly rich example of a non-Markovian process in 
which the bias at each site disappears when all the cook- 
ies on the site are eaten. Because a random walk in one 
dimension is recurrent, that is, the probability to even- 
tually return to the starting point equals one, the walk 
returns to previously-visited sites infinitely often. For 
the 1-excited walk, the reinforcement mechanism due to 
the presence of cookies is insufficient to counterbalance 
this recurrence of the random walk; thus the 1-excited 
random walk remains recurrent for any p < 1. Never- 
theless, this model exhibits many novel features. For 
example, the probability that a 1-excited walk on an in- 
finite line first returns to the origin at time t decays as 
^ '■^"^ rather than the t~^^^ decay for a pure random 
walk UllSIIil- 

On the other hand, the fc-excited random walk with 
fc > 1 is more strongly influenced by the cookie-induced 
bias. In fact, for p > ifl -I- i), the fc-excited random 
walk becomes transient ]a| • This transition has dramatic 
consequences on the statistical properties of the walk as 
p is varied. As a final introductory point, the excited ran- 
dom walk appears to be most interesting in one dimen- 
sion because of the non-trivial competition between the 
reinforcement when a cookie is eaten and the fact that 
encounters with sites that contain cookies become pro- 
gressively less frequent. In contrast, an unbiased random 
walk in greater than one dimension typically encounters 
previously- unvisited sites Because each such site 

contains a cookie, the reinforcement mechanism acts as 
a non-zero global bias that is superimposed on the ran- 
dom walk motion. Thus the asymptotic behavior of an 
excited walk in higher dimensions is close to that of a 
simple biased random walk. 

In the next section, we begin by recounting some ba- 
sic first-passage properties for diffusion in a finite inter- 
val . We present these results separately because 
they underlie our analytical treatment of excited random 
walks. In Sec. Ill, we study a particularly simple version 
of the excited random walk - the 1-excited walk with a 
refiecting boundary condition at the origin. We use ba- 
sic tools of first-passage processes to determine the time 
dependence of the size of the region in which the cookies 
have been eaten and thence the probability for the walk 
to first return to the origin. The long-time behavior of 
the first-passage probability in a closely related model 
was obtained previously, using quite different methods, 
in 0, Kill . Our approach is also complementary to the 
tools employed in the probabilistic studies of excited ran- 
dom walks given in Ref. naiii. In Sec. IV, we first 

present numerical simulations for the unusual properties 
for the probability distribution and for the span of the 
1-excited random walk on the infinite line. We also con- 
struct a probabilistic argument to quantify the asymme- 
try in the number of cookies eaten from the left and right 
sides of the line. Finally, in Sec. V we present several in- 



triguing properties of the 2-excited walk on the infinite 
line, including the lack of scaling in the probability dis- 
tribution, the time for the 2-excited walk to traverse a 
finite domain, and basic features of the cookie density 
profile that is left behind as the walk advances. 



II. PRELIMINARIES: FIRST PASSAGE IN AN 
INTERVAL 

To determine the motion of excited random walks, we 
require the first-passage properties of a nearest-neighbor 
unbiased random walk, or equivalently, continuum diffu- 
sion, in a finite interval [0, k] when both boundaries are 
absorbing. The first such property is the eventual hitting, 
or exit, probability to a; = and to a; = fc for a random 
walk that starts at xq within the interval [Tl. IT^ . These 
exit probabilities are: 



E.{k,xo) = l-^ 



(1 \ ^0 



(1) 



Closely related quantities are the conditional exit times 
to X = and x — k. These are the mean times for 
the random walk to hit a; = or to hit x = k with the 
restriction that the other boundary is never hit. These 
conditional times are: 



t^(fc,xo) = ^(2fc^ 



-xq) t+{k,xo) = ^{k'-xl) . (2) 



Next, we will also need the Laplace transforms of the 
first-passage probabilities for a random walk that starts 
at xo to exit at the left edge and at the right edge of 
the absorbing interval [0, fc] at time t. In the continuum 
limit, these first passage probabilities are (Eq. (2.2.10) in 
Ref. Ill): 



£k{xo, s) 

Tkixo, S) 



sinh \/2s {k — xq) 



sinh v2s fc 
sinh \/2s xq 



sinh V 2s fc 



(3) 



where we take the diffusion coefficient to be D = 1/2 for 
the nearest-neighbor random walk. 

We will also need first-passage properties in the inter- 
val [0, fc] when the boundary at x = is refiecting and 
the other boundary is absorbing. We implement the re- 
fiecting boundary condition for a discrete random walk 
that is at site x = by allowing the walk to remain at 
X = with probability 1/2 or by having hop to a; = 1 
with probability 1/2. This definition for the reflecting 
boundary has the advantage that it leads to a constant 
probability density around this boundary for large times. 
With this deflnition, the average time to exit at the right 
edge of the interval [0, fc] as a function of xq is 



t(fc, Xo) ~ k{k + 1) — a;o(a;o + 1) 



(4) 



Finally, the Laplace transform of the flrst-passage proba- 
bility to exit at the right edge at time t in the continuum 
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limit is (Eq. (2.2.21) in "l^ 



7?.fc(xo,s) 



cosh v2s (k — xq) 
cosh \/2s k 



(5) 



We will exploit these results to now determine various 
basic features of the 1-excited random walk in one di- 
mension. There are two basic properties of the excited 
random walk that we will be able to determine analyt- 
ically: (i) the location of the last (leftmost) cookie; we 
obtain this property in terms of the first-passage statis- 
tics of a 1-excited random walk to the right edge of the 
cookie-free domain, and (ii) the statistics of returns to 
the origin. 



III. 1-EXCITED RANDOM WALK ON THE 
HALF-LINE 

Consider the 1-excited random walk that is restricted 
to the non-negative integers by imposing a reflecting 
boundary condition at the origin. Initially there is one 
cookie at each site and the walk is at n = 0. At the 
first time step, the walk either hops to n = 1 with prob- 
ability p or it remains at site n = with probability 
q (Fig. At later times, a bias in the random walk 
motion occurs whenever the walk reaches the right edge 
of the cookie-free domain. The restriction of fixing one 
boundary condition while the other boundary condition 
"floats" greatly simplifies the analytical description of the 
1-excited walk on the half-line compared to the infinite 
line system. 



A. Time to Hit the Last Cookie 

It is useful to characterize the walk by the sequence of 
times tn at which the walk hits the leftmost (last) cookie 
at n (Fig.^. Suppose that the walk has just reached the 
cookie at n — 1. Then with probability p, the next step 
of the walk is to right, which takes one time step, while 
with probability q, the walk steps left to n— 2 and the last 
cookie now is at n. In this latter case, the average time it 
takes the walk to reach site n is t{n,n — 2) = An — 2 [see 
Eq. plus the time for the initial step from n — 1 to 
n — 2. Thus the mean time increment for the last cookie 
to move from n — 1 to n is 

{dtn) = p-l + q-{l + 4n-2) = l + q{4n - 2) . 

Consequently, the mean time for the walk to hit the last 
cookie at n for the first time is 



(tn) =^(rftfc) =2qn^ 



(6) 



k=l 



Notice that for g ^ (p — > 1), the early time behav- 
ior is governed by ballistic motion, after which there is a 



crossover at a time of order 1 /q to the asymptotic diffu- 
sive behavior of {tn) ~ 2gn^. 

We now extend this analysis to obtain the distribution 
of first-passage times for the walk to hit the last cookie 
at n at time t. Let Qk{t) be the probability that the 
time interval between the walk first hitting site fc — 1 and 
first hitting site k equals t. Similarly, let Fn{t) be the 
first-passage probability that the walk hits site n for the 
first time at time t. Then 



Y[Qkitk)dtk. 



k=l 



with ti + t2 + + ■ ■ ■ tn = t. Because this relation is a 
convolution, the Laplace transform Fn{s) obeys 



Fnis) = l[Qkis). 



(7) 



k=l 



Now the Laplace transform of Qkit) is, using Eq. (O, 

Qk{s) = p(l-s)+(?(l-s)7^fc(2,s). (8) 

The first term accounts for the direct step from A: — 1 to 
k. The contribution of this direct process to the discrete 
generating function would be pz. In the continuum limit, 
the factor z maps to the Laplace variable 1 — s. The 
second term accounts for a single step to the left (factor 
q(l — s)) and the subsequent first-passage from A: — 2 to 
k in the presence of a reflecting boundary at the origin. 
The first step of the walk is accounted for by defining 
Qi{s) = e'. 

We anticipate that the dominant contribution to Fn{s) 
in the long-time limit comes from terms where k ^ \Jt 
[T^. In this hmit, Eq. ((H)) reduces to 

Qfe(s) = 1 - 2q\/2stanh\/2sfc , 

by expanding 7?.fc of Eq. jSJ for large fc. Using this ap- 
proximation for Qk and taking the limit s — > 0, corre- 
sponding to long-time behavior, Fn{s) can be written as 

n 

Fn{s) « Jl (l - 2gV2itanh\/2iA;^ , 

fc=i 

= exp I ^ ln(l - 2qV2s tanh V2s fc) J ' 



(9) 



exp [ — / 2gv2stanhv2sfc dk 
(cosh\/2ia;)"29 ^ 



where we use / tanh xdx — In cosh x and we replace the 
discrete variable n by the continuous variable x in the 
last two lines. 

From this formal solution, the cumulants Kg of the first- 
passage time may be obtained as the coefficients of (— s)^ 
in the power series of ln_F(x, s). For this we use the 
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product representation of the cosh function, and the se- 
ries expansion of ln(l + z) to give 

lnF(x, .) = 2,^ i ( ^ ) i-sf Y^{2k + 1)-^ . 
i=i \ / k=0 

Recognizing that the second sum can be expressed in 
terms of the Riemann zeta function, the cumulants are 



(10) 



where 



^)(£-l)!7r-^^C(2£), 



which gives C(l) = 2,C(2) = 4/3, C(3) = 32/15, etc. 
This resuh agrees with the direct solution given in Eq. © 
for the average first-passage time {t{x)) = ki — 2qx^. 

Finally, the time dependence of the first-passage prob- 
abihty for the last cookie - the probability that the cookie 
at X is first hit at time t - is obtained by inverting the 
Laplace transform in Eq. I^. To simplify this calcula- 
tion, we introduce the scaling variables 



X 

VI 



St , 



to rewrite the first-passage probability F{x,t) in the scal- 
ing form 

tF{x,t) = <f{pi) J dy e^(coshy2^A*)"'^ • (H) 



For the case of no bias {p = q = 1/2), the integrand 
has simple poles on the negative real axis at 



y 



(2/c + l)27r2 



fc^O, 1,. 



By evaluating the corresponding residues at these poles, 
we find the scaling function 



k=0 



(12) 



Because of the presence of an exponential cutoff, this 
series is easy to compute numerically and the result is 
displayed in Fig. El together with the small-/^ limiting 
behavior (corresponding to small x for fixed t or large t 
for fixed x) that stems from the first term of Eq. H12|l : 
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FIG. 2; The scaling function ^(/i) in the scaling limit for 
p = 0, 0.5, and 0.9. The large-/! asymptotic is represented by 
the first term in Eq. 1141 (thin curves). For p — 0.5 the first 
term of the small fi series of Eq. I12II is also shown (dotted) . 



Each term in this series can be easily back transformed 
term by term (see 29.3.82 in Ref. 16]). The scaling func- 
tion thus has the form 



*(m) - 



22<;+l/2 



■E 

fc=0 



-2q 
k 



{k + q)e 



(14) 



For q = 1/2, this function is identical to <i>(/i) in Eq. H12|l 
in spite of the apparent visual differences in these two 
formulae, since both functions originate from the same 
Laplace transform. Eq. H14|l can be viewed as the large- 
/i expansion of $(/x). As shown in Fig.|21 the k = term 
in this series representation 



« ^lq 



22g+l/2 



(15) 



together with the fc = term in the small-^ series for 
Eq. H12|l give a fairly good approximation for $(/x) over 
the entire range of for p = 1/2. Fig. |3| shows the so- 
lution Eq. H14() . together with the exact numerical prop- 
agation of the probability distribution at different times 
(see Appendix). 



B. Location of the Last Cookie 



For general p ^ 1/2 the poles in Eq. (|ll|l are branch 
points and we therefore evaluate this integral in a sim- 
pler manner by first using the binomial theorem to write 
Eq. jnj as 



fc=0 



-.(k+q) 



(13) 



Another important characterization of the excited ran- 
dom walk is the probability distribution for the location 
of the last cookie. Now the last cookie either remains 
stationary or it moves only to the right during each step 
of the underlying excited random walk. Consequently, 
we can write a simple relation between the first-passage 
probability for hitting the last cookie, derived in the pre- 
vious section, and the probability distribution of this last 
cookie. 
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FIG. 3; The scaling function compared to the numerical exact 
propagation of the probability distribution. 



The probability that the last cookie is located at site 
n at time t equals the joint probability that a walk first 
hopped to site n — 1 before time t, and that this walk 
then first hopped to site n after time t. That is, 

t oo 

Pn{t)^Y. E Fn-lit')Qn{t") . (16) 
t' = l t"=t-t' + l 

In the continuum limit, the Laplace transform of this 
equation is 

Puis) = -F„_i(s) [1 - Q„(s)] =. - [F„^i{s) - F„(s)] , 
s s 

where we used Eq. ((7J to simplify the last factor. In the 
continuous space limit, the relation between P and F 
becomes 



P{x,s) = - 



1 dF{x,s) 
s dx 



which is equivalent to 



P{x,t) = ~-^ f dt'F{x,t') ■ (17) 
ox Jq 

Substituting the scaling form for F{x,t) into Eq. H17|) 
and performing the derivative, we express the integrand 
in scaled units, after which the integral can be evaluated 
easily. 

These steps lead to the following simple connection 
between the first-passage and occupation probabilities of 
the last cookie: 

F(x,i) = 7*(4), and P(a:, i) = -$(^) , (18) 

with the same scaling function for both quantities. We 
exploit this connection between F and P to obtain the 



finite time data in Fig. by first calculating the position 
distribution of the last cookie from the numerical evalu- 
ation of the Master equation (see Appendix), and then 
using Eq. (|18|l to arrive to the scaling function $(/i). 

Finally, we use Eqs. H14|l and H18() to obtain the exact 
expression for the average position of the last cookie: 

/•oc 

{x{t)) ^ dx xP{x,t) = 2VtAf(q) , (19) 
where the proportionality constant 

Jo V^Q r(H-2g) 

is obtained by first integrating (|14|l term by term and 
then using the identity 



E 

k=Q 



2q\ 1 [r(i + <?)]^ 



k J k + q qT{l + 2q) 



The surprising point of this calculation is that {x{t)) ^ 
Vt/q as g — *■ 0, while the mean first-passage time for the 
walk to reach a distance x scales as {t{x)) cx qx"^. More 
precisely. 



{t{x)) = 2qx^ 



(xit)) 



1 



/2Trq 



Vt , (20) 



while naive dimensional scaling would suggest that (t) cx 
qx^ implies {x{t)) cx y/tjq. The fact that this naive 
scaling does not hold is a result of the broadness of the 
distribution of first-passage times for the last cookie. 



C. First Passage to the Origin 

As a complement to the statistical properties of the last 
cookie, we now study first passage to the origin, namely, 
the probability for the 1-excited walk to reach x = 
for the first time at time t. We first provide a heuristic 
argument (previously given in Ref. 0) that the 1-excited 
walk is recurrent and then we extend this line of reasoning 
to argue that the first-passage probability asymptotically 
scales at t~'i~^. 

If the walk is at a point x — 1 > for the first time, 
then the probability for the walk to subsequently reach 
X without hitting the origin \s p -\- q{l — = 1 — 
Here the factor p accounts for the event where the walk 
immediately eats a cookie by making a step to the right, 
and the factor 9 (l — f ) [see Eq. ^] accounts for the 
probability that the first step is to the left and that the 
walk eventually eats another cookie on the right without 
hitting the origin. Thus, if the walk starts at site x ~ 1 
with the system initially full of cookies, the probability 
for the walk to reach (at least) an arbitrary point x > 1 
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without hitting the origin is 



P{x) 



n(l-|)=exp 

" 2q 



k=2 



.k=2 



exp 



■dy 



(21) 



Since a;"^"? — > as a: — *■ c», the random walk must neces- 
sarily return to the origin; that is, the 1-excited random 
walk is recurrent, unless p — \. 

Since P{x) is the probability that the walk eats at least 
X cookies without hitting the origin, P'{x) ^ x~-^~'^'^ is 
the probability that the walk eats exactly k cookies. The 
average number cookies eaten by a walk that does not 
hit the origin is therefore (n) ^ J dx x~^'^. Thus for 
q > 1/2 (p < 1/2) a walk that manages to avoid the 
origin typically wanders only a finite distance away (for 
a more rigorous argument see Sec. lIV(j)l . 

If the walk is at site x — 1 for the first time, the con- 
ditional average time to reach site x without hitting the 
origin is one with probability p, and |(a; — 1) with prob- 
ability q (see Eq. (0)), that is {dtx) = I + ^{x - 1). 
Thus the average time to reach site x, without hitting 
the origin, is (t^) = 1 + J^t^si^^k) ^ for large 
x. From this relation between t and x, we infer that 
the survival probability, namely, the probability that the 
walk does not reach the origin by time t, should vary 
as Sit) = P{x) ~ x-^? - Consequently, the first- Then 



passage probabihty F{t) 
decay as 



as 

dt 



should asymptotically 



Fit) ^ t 



(22) 



For q < I, the average time to reach the origin is infinite. 
Quite strangely, although it takes infinite time for the 
walk to hit the origin for any q, the walk eats only a finite 
number of cookies, i.e., wanders only a finite distance 
from the origin, for q > 1/2. This result remains valid 
even if there is an initial cookie-free region at the origin, 
as this initial condition would affect the pre-asymptotic 
factors in Eq. (ED- 

To compute the first-passage probability to the origin 
more precisely, it is useful to define two auxiliary quan- 
tities, Lk{t) and Rk{t), as follows: 

• Lfe(i) =the probability of first hitting fc = without 
hitting fc + 1 a time t after the walk first hits k. 

• Rk{t) = the probability of first hitting fc -I- 1 without 
hitting fc = a time t after the walk first hits k. 

In the Laplace domain, these first-passage probabilities 
are, using Eq. Q, 

Lk = (7(1 - s)4+i(fc - 1, s) , 

Rk = p[l- s) + q{l- s)rk+i{k -l,s) . 



For example, for L^, the first step must be to the left 
(factor q{l — s)) and then the factor ik+i is the first- 
passage probability to the origin, without hitting fc + 1, 
for a walk that starts at fc — 1. Similarly, for Rk, the term 
p{\ — s) accounts for the transition fc — > fc -I- 1 in a single 
step. The factor q{l — s) in the second term accounts 
for a single step to the left after which the walk is at 
fc — 1 while the last cookie is at fc -I- 1. Then the factor 
rfc+i(fc — 1, s) gives the first-passage probability to fc + 1 
without hitting the origin, when starting from fc — 1. 

Finally, we obtain the Laplace transform of the first- 
passage probability to the origin for a walk that begins at 
X = 1 with the system initially full of cookies by summing 
over all paths that contain 0, 1, 2, . . . first passages to the 
last cookie before the origin is reached. This gives 



F{s) = Li+Ri[L2 + R2{L^ + R^{U + ... 

= il + R1L2 + R1R2L3 + R1R2R3L4 



(23) 



In the limit s ~ 0, this Laplace transform is just the 
integral of the first-passage probability over all time or, 
equivalently, the probability of eventual return to the ori- 
gin. It is easy to verify that this return probability equals 
one. In the s limit, the auxiliary probabilities Rk 
and Lk are 



_2q_ 

k + 1' 



_2q_ 

k+l 



1-Rk. 
(24) 



F{s = 0) = 1 - Ri + Ri{l - R2) + RiR2{l - R3) 

+ RiR2R3il - R4) + ■ . . . (25) 

This expression obviously equals one for all p < 1, while 
for extreme case of p = 1, F{s = 0) = 0. Thus eventual 
return to the origin is certain for any bias p < I. 

More formally, we denote the sum of the first m terms 
of Eq. as 



m k—1 

fc=i j=i 



(26) 



For finite m, we use Eq. (|24|l to give, in the limit s — s- 0, 
2q 



Fmi-S = 0) 



r(2 - 2,) 1^ r(* + 2 



which goes to one as m ^ 00 . The use of the s = 
forms for Rk and Lk from (|24l) remains valid for m of 
the order of or less. We now assume that the main 

contribution to F(s) comes from terms with m < l/^/s 
in Eq. (jSEJ- That is, 



F{s) 



2q 



r(2-2g)^ r(fc-t-2) 
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This finite sum can be written as the difference of the 
infinite sum (which equals 1) and the sum from l/\/s to 
oo. In this second sum, fc ^ 1 for all terms as s oo. 
Therefore we can replace the ratio of gamma functions 
by k°'~'' and the sum by an integral. These steps lead to 



l~F{s) 



dk k 



-2q~l 



(29) 



This result for F{s) implies that for large times the 
first-passage probability to the origin indeed behaves as 
Fit) ^ t- 



-l-q 



IV. 1-EXCITED WALK ON THE INFINITE 
LINE 

A. First-Passage Probabilities 

Using the same first-passage analysis as that given for 
the 1-excited walk on the half-line, we can immediately 
deduce the first-passage probabilities for the 1-excited 
walk on the infinite line to hit the left side and the right 
side of a cookie-free region. Thus the first-passage prob- 
ability for the walk to hit the left side of this gap is the 
same as the first-passage probability to the origin for the 
half-line system. From Eq F^it) ^ where 

we now write the subscript L to denote hitting the left 
side of the gap. The first-passage probability to the right 
side of the gap can immediately be obtained merely by 
interchanging p and thus Fji{t) ~ t^^^'P. Since the 
exponents in both of these first-passage probabilities are 
smaller than 2, the mean time to hit either edge of the 
cookie-free domain is infinite. 



B. The Probability Distribution 

Another basic characteristic of the 1-excited walk is the 
probability distribution of displacements. An analytical 
treatment of this quantity is more involved than the semi- 
infinite system because of the existence of two moving 
boundaries. We thus turn to numerical simulations to 
provide a qualitative picture of this distribution. 

Fig. 0] shows the scaled probability distribution for 
p = 0.7, 0.8, and 0.9. For p > 0.5, the probability distri- 
bution has a nearly Gaussian, but slightly asymmetric, 
appearance. As p increases toward 1, this asymmetry 
becomes more pronounced, as shown in the figure. In 
addition, an apparent singularity in the distribution de- 
velops near the origin as p — > 1. By running simulations 
for p = 0.7 and 0.95 to 431, 440 « 1.5^^ steps, we observe 
quite accurate data collapse for p = 0.7. However for 
p = 0.95, there is a small but persistent departure from 
data collapse for < x/a < 0.3 (inset to Fig. ^ We 
do not understand whether this anomaly arises from a 
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FIG. 4: Scaled probability distributions of the 1-excited ran- 
dom walk for 2 x 10® realizations for 85,223 « 1.5^* time 
steps for p — 0.7, 0.8 and 0.95. Plotted on the horizontal axis 
is x/a and on the vertical axis is aP{x,t), where P{x,t) is 
the probability distribution of the walk. The data have been 
averaged over a ^% range. The inset shows the distribution 
for p — 0.95 near the origin for still longer times: t = 1.5*, 
with fc = 29 - 32 (127,834 - 431,440 steps). 



slow correction to scaling term in the probability distri- 
bution or whether there exists an effect that lies outside 
of scaling. 

For all values of p > 1/2, both {x{t)) and the variance 
a{t) = \/{x{tY) - {x{t)Y grow as t^/"^. The time de- 
pendence for the displacement is easy to understand on 
heuristic grounds. Once a cookie- free domain of linear di- 
mension L has been created, the domain then grows only 
when the walk reaches either edge of the domain. Since 
the motion of the walk is purely random within the do- 
main, the probability that the walk hits either edge is 
proportional to 1/L. If the walk hits the right edge of 
the domain, then the walk will move ballistically for typ- 
ically 1/(1— p) = 1/q steps as p ^ 1. Thus the equation 
of motion for the domain length will be i = 1/ [Lq) as 
p ^ \. The solution to this equation gives L cx t^^"^ for 
any value of p < 1. However, it bears mention that this 
naive approach gives the wrong value for the amplitude 
- 1/ Jq - rather than 1 /q that was derived in Eq. (|20(l . 



C. Asymmetry 

A useful geometric characteristic of the 1-excited walk 
is the asymmetry in the span of the walk. Here the span 
is the length of the cookie-free region. For an unbiased 
random walk, the growth of the span with time and its 
distribution are well understood [ig, and we now inves- 
tigate how the bias inherent in the 1-excited walk influ- 
ences the span. It is useful to characterize the span by 
the positions of the last cookie on the right and the last 
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cookie on the left of the cookie-free region. We can un- 
derstand the motion of these extreme cookies by the fol- 
lowing heuristic argument. Consider the situation where 
the cookie-free gap has length L — 2 and where the walk 
has just stepped to the extreme cookie on the right. As 
shown at the beginning of Sec. lIIld the probability that 
the walk eats another right cookie without eating a left 
cookie is 1 — 

With this basic result, the probability that the walk 
eats precisely r > consecutive cookies (we term this 
event a single "meal" ) from the right edge of the cookie- 
free region is 



r-2 



p{r) = n 1 



_2q_ 

L + J 



2q 



L + r-1 



j=o ^ (30) 
r(L) T{L + r-l-2q) 



2q 



T{L-2q) T{L + r) 



where we use the convention Yl'^a ^ ^ ^'-"^ a > b. The 
probability that the walk eats precisely one cookie (the 
one at the starting position) is then 2q/L. 

When the initial cookie-free gap is large, this probabil- 
ity simplifies to 

P{r) r^2qL'^i {L + r)-'^i-^ . 

The L dependence can be scaled out using the dimen- 
sionless variable f ^ r/L and the rescaled probability 



P(f) = LP{r) = 2q{l + r)-^'^-^ 



(31) 



Note that this probability is normalized, drP{r)=l. 
Using this probability, the average relative number of 
consecutive cookies eaten from the right side of the gap, 
Jp°° df rP{f), is 

{f} = {r/L) = \ 2q-V (32) 
I oo, q<l/2. 

This result is paradoxical at first sight because, as we 
showed in Sec. IIII CI the average time to return to the 
left side of the gap is infinite. However, the walk eats, on 
average, only a finite number of cookies if q < 1/2, even 
though the time taken for its meal on the right side of 
the gap is infinite. 

The average relative number of eaten cookies from the 
left (/) is also given by Eq. H32|) after interchanging q 
and p. Although the 1-excited walk eats, on average, 
infinitely many consecutive cookies from a single side of 
the gap during a single meal, the average ratio {r/{l + r)) 
turns out to be finite for any value of q. To show this 
we compute the joint probability of eating r consecutive 
cookies on the right and then I consecutive cookies on 
the left before the next right cookie is eaten, under the 
initial condition that the walk starts at the position of 



the right extreme cookie. This joint probability can be 
derived analogously to Eq. Ip!n| 

T{L + I - 1 - 2p)T{L + I + r - 1 - 2q) (33) 

V{L + l-2q) T{L + l + r) 

- 4g (1 - q)L^''{L + rf-^'i{L + r + if'^. 

It is again useful to define the rescaled variables I = l/L 
and f — r/L, as well as the rescaled probability 

Pir, S) = T^P{r, I) = Aq{\ - q){l + ff-^^l + Sf^-^ , 

(34) 

where s = I -\- f is the total number of cookies eaten in 
one right-left cycle. One may check that this probability 
is normalized, ds drP{r, s) = 1. 




FIG. 5: Ratio of the number of right cookies eaten to total 
cookies eaten after a sine le cycle (EHJ- Tlie squares are ob- 
tained from simulations with an initial gap ranging from 20 
to 200. 

Finally, the average ratio for the number of right cook- 
ies to total cookies eaten in a single cycle is 



ds I dr — P{r, s) 
Jo ^ 



^Aq{l-q) ^il + sf^' 

Jo s Jo 



dr r(l+r) 



l-4q 



(35) 



This integral can be performed with the help of 3.235 in 
Ref. |l?| (both for g < 3/4 and g > 3/4), and also using 
some basic properties of the Euler ip function J^] to give 



^\ = l--^--^^^%^cot2.g (36) 
s/ (l-2g)2 (l-2g)(3-4g) ^ ^ ' 



This gives the asymptotically exact ratio for the number 
of right cookies to total cookies eaten in one cycle, under 
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the assumption that the initial cookie-free gap is large. 
Although (l) or (r) is infinite, the g-dependence of the 
ratio {r/{l + r)) is well-behaved (Fig.EJ. 

Despite the fact that the average number of consecu- 
tive cookies eaten from the left or the right side of the 
gap in one cycle is infinite, our data clearly show that the 
positions of the two extreme cookies both grow as t^/^, 
but with very different amplitudes. To obtain an approx- 
imation for these amplitudes we now calculate the typical 
number I* and r* of eaten consecutive cookies from each 
side. 

We estimate the value of r* by setting the probability 
of eating at least r* consecutive cookies on the right equal 
the probability of going to the left before r* right cookies 
are eaten. That is. 



n 



1 - 



2q_ 

L + J 



(37) 



This criterion gives a crude estimate for the number of 
right cookies eaten during a single consecutive string. For 
L » 1, Eq. ^\ becomes 



L 



2q 



1 

2 ' 



which leads to the estimate that the typical length of the 
string of eaten right cookies is r* — (2^/^'' — 1)L. Ap- 
plying this same argument to the left side of the cookie- 
free region, a string of I* — (2^/^^ — 1)L consecutive left 
cookies are typically eaten. Their ratio thus equals the 
constant value 



21/29 _ 1 
2l/2p _ i' 



(38) 



As a function of time, r* /I* clearly approaches a constant 
value as t — > oo and it is easy to estimate its asymptotic 
value which is shown in Fig. El together with the predic- 
tion of Eq. (PI). 

To conclude this section, we point out another intrigu- 
ing feature of the extreme cookies; namely the probabil- 
ity distribution in their positions. Quite strikingly, these 
distributions have short-distance power-law tails, both 
for the left and for the right extreme cookies (Fig. ^ 
for p 7^ -i. In contrast, for the pure random walk case 
of p = i, the extreme cookie distribution is essentially 
the same as the span distribution, which is known to be 
a Gaussian- like function The origin of the short- 

distance power-law tails in the probability distributions 
of the last cookie seem to be connected with the anoma- 
lous scaling of the first-passage probability, but we do not 
have a definitive explanation for this phenomenon. 



V. 2-EXCITED RANDOM WALK 

In the 2-excited random walk, each lattice site initially 
has two cookies and the (p, q) bias occurs as long as the 




FIG. 6: Tlie ratio r* /l* as a function of p for the 1-excited 
random walk. The smooth curve is the the prediction from 
Eq. 13811 . while the data are based on 2 x 10^ realizations of 
the walk. 
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FIG. 7: Probability distributions for the position of the right- 
most cookie (positive slopes for small x) and the leftmost 
cookie for the 1-excited random walk. Data are based on 
2 X 10® realizations for 85, 223 time steps, for p — 0.7 (line), 
0.8 (dashed), and 0.9 (dotted). For the rightmost cookie, 
the apparent power law for small x has exponent 0.35, 0.54, 
0.74 and 0.84 for p=0.7, 0.8, 0.9 and 0.95, respectively, while 
for the leftmost cookie, the respective exponents are —0.42, 
—0.63, —0.83, and —0.94. Also shown is the data for the 
unbiased random walk (o). 



walk lands on a site that contains at least one cookie. 
This 2-excited random walk undergoes a transition be- 
tween recurrence and transience as p is increases beyond 
a critical value Pc = | @ ■ Indeed, the fc-excited walk for 
any k > 1 appears to have qualitatively similar behavior 
as the 2-excited walk. For simplicity, we thus focus on 
the 2-excited random walk. 

We first recount the qualitative argument for the tran- 
sition between recurrence and transience of the fc-excited 
walk that was given in Ref. ^ . Whenever the walk is in 
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a region that contains cookies, there is ostensibly a bias 
velocity v = p — q. Thus the time that the walk takes to 
travel a distance L should he L/{p — q). If this time is less 
than the total number of cookies kL in this region, then 
the walk does not have sufficient time to consume all the 
cookies as it moves. Consequently, the bias will persist 
forever and the walk will be transient. The condition for 
transience is therefore L/v < kL, or 

P>P. = i(l + ^). (39) 

As we shall see, the behavior of the 2-excited walk is 
qualitatively different in the regimes p < Pc and p > Pc- 



1.2 




x/a(t) 



A. Probability distributions 

Fig. |S1 shows scaled displacement probability distribu- 
tions for the 2-excited random walk for representative 
values of p. As p passes through pc, there is a sud- 
den drop in the integrated probability that the walk 
lies in the range x < 0. Another basic feature of the 
distribution is the departure from scaling, as shown in 
Fig. 1^1 For p = 0.9, data collapse has not yet set in by 
431,400 « 1.5'^^ time steps. If scaling were to hold, it 
would only arise after an extraordinarily long time. In 
the absence of such a clearly identifiable long time scale, 
it seems that the asymptotic distribution will not obey 
single-parameter scaling. 




x/a(t) 

FIG. 8: Scaled probability distributions of the 2-excited ran- 
dom walk for 85,223 time steps and for 2 x 10® realizations. 
The data have been averaged over a i% range. 



FIG. 9: Scaled probability distributions of the 2-excited ran- 
dom walk for p = 0.9 at different times for 2 x 10® realizations. 
The data have been averaged over a ^% range. 

is shown. For p — 0.6 this density profile follows single- 
parameter scaling and the also decays to zero near the 
origin because of the recurrence of the walk for p < 3/4. 
On the other hand, for p = 0.9, the cookie density shows 
a very small departure from scaling and the density goes 
to a non-zero limit near the origin as a result of the tran- 
sient nature of the walk. We can understand this limiting 
cookie density near the origin by the same argument that 
led to the criterion for pc in Eq. H39|l . For general p > Pc, 
a walk typically resides for a time l/{p — q) at each site 
and should cat l/{p— q) cookies, thus leaving a residue 
of 2 — l/(p — g) uneaten cookies per site. For the case 
p = 0.9 shown in Fig. 1101 we therefore expect 3/4 of a 
cookie per site near the origin. 
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The transition in the behavior of the probability dis- 
tribution is mirrored in the spatial profile of the uneaten 
cookies. For p < pc, the 2-excited random walk is recur- 
rent and thus the region close to the origin should ulti- 
mate be devoid of cookies. This behavior is illustrated 
in Fig. ^1 where the scaled density of uneaten cookies 



FIG. 10: Cookie profile for the 2-excited random walk at p = 
0.6 (lower data) and at p = 0.9 (upper data) for times 656, 
1477, 3325, and 7480. The x-axis has been scaled by the 
average displacement of the walk at time t. 

A more detailed characterization of the cookie profile 
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is provided by the size distributions of domains that con- 
tain one cookie and contain zero cookies. The distribu- 
tion of one-cookie domains is an exponentially decaying 
function of domain length (Fig. Ill|l . This behavior is 
simple to understand. The process that gives rise to 
a one-cookie domain as p 1 is a walk that makes a 
single step to the left after a run of consecutive steps 
to the right. During this run of rightward steps into 
previously-unvisited territory (each occurring with prob- 
ability p), one cookie is left behind at each site. A step 
to the left (which occurs with probability q) terminates 
this one-cookie domain by a cookie-free gap. Thus the 
probability of a one-cookie domain of length k is simply 
p'^q — ge^'^''"(^/P^. From linear least-squares fits to the 
data in Fig. EI the slopes for p = 0.6, 0.7, 0.8, and 0.9 
are, respectively, -0.472, -0.348, -0.219 and -0.104. 
On the other hand the corresponding slopes from by the 
distribution (J e-'^'^^^i/p) are -0.511, -0.357, -0.223 and 
—0.105, in excellent agreement with the data, especially 
as p — > 1. 




20 40 60 80 100 120 140 
length 



FIG. 11: Length distribution of one-cookie domains on a semi- 
logarithmic scale for 50000 realizations at time t = 37, 877 « 
1.5^6 for p = 0.6, 0.7, 0.8 and 0.9. 

The length distribution of zero-cookie domains is more 
intriguing. Such domains arise both by the growth of 
isolated zero-cookie domain and by the merging of two 
neighboring zero-cookie domains when an intervening 
one-cookie domain disappears. The zero-cookie domains 
also have qualitatively different behavior for p < pc and 
p > Pc- In the former case, the distribution has a sec- 
ondary peak close to the maximum- length domains, while 
the distribution appears to have a pure power-law form 
for p > Pc. 



B. Transit time 

Our numerical simulations indicate that the mean dis- 
placement {x{t)) of the 2-excited walk grows with time 
as , with v apparently dependent on p. For the 3 cases 




length 

FIG. 12: Length distribution of zero-cookie domains on a 
double logarithmic scale for the same parameters as Fig. 1111 

of p = 0.7, 0.8 and 0.9, least-squares fits to the data in 
the range t > 20000 give i/ = 0.54, 0.64, and 0.85, respec- 
tively. Naively, one might anticipate that the transience 
of the walk for p > 3/4 would lead to an average dis- 
placement that grows linearly with time. Indeed, the ar- 
gument that gave the threshold bias = 3/4 was based 
on the walk possessing a non-zero bias velocity propor- 
tional to p — q for p > Pc- On the other hand, it has been 
proven that the mean velocity of the 2-excited walk, de- 
fined as \imt^oc{x{t)) /t equals zero for all p < 1 0. 

While we do not know how to determine the time de- 
pendence of the displacement, we can address the com- 
plementary problem of estimating how long it takes a 
2-excited walk to traverse a finite interval. We already 
know, from the 1-excited walk, that the mean time to 
traverse any cookie-free interval within an infinite sys- 
tem is divergent. Thus we now consider the time for a 
2-excited walk to traverse a finite cookie-free interval. As 
we shall show, the transit time grows faster than linearly 
in the interval length when p is close to but strictly less 
than 1. 

The 2-excited random walk starts at x = and we 
compute the time to reach x = L, where {1 — p)~''. 
This inequality ensures that a "defect" , namely, a region 
where the walk takes a small number k of consecutive 
steps to the left after a long string of rightward steps, will 
certainly occur in the interval. This defect represents the 
lowest-order perturbation to the purely ballistic motion 
for the 2-excited walk when p ^ 1. After the defect has 
been completed and a step to the right has eventually 
been made, the walk finds itself at the point x = 2 in a 
cookie-free region [0,k]. We ask: what is the mean time 
for the walk to reach the right edge of this interval? As 
illustrated in Fig. El this system is a 1-excited random 
walk in disguise. 

We denote by the mean time for the 2-excited walk 
to reach the point x = k in this geometry. To simplify 
the calculation, we assume that whenever the left edge 
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FIG. 13: A 2-excited random walk in which the walk has made 
k consecutive steps to the left after a long run of rightward 
steps. The corresponding cookie configuration is also shown. 

of the domain is reached and the last cookie is eaten, the 
next step of the walk is always to the right. Because this 
next step could be to the left, our approach represents a 

lower bound to the transit time that we denote by T^: 

As defined in Eq. and lO, with probability 
-B-i-(/c, 2), the walk reaches x = k without touching a; = 
and the time for this event is t-|- (fc, 2) . On the other hand, 
with probability E^{k, 2) the walk hits the left edge and 
then, according to our approximation, hops one step to 
the right. Since this second argument always equals 2 
in all formulae involving E± and t±, we henceforth drop 
this argument for simplicity. After these two events, the 
walk is a distance 2 from the left edge of a cookie-free 
domain of length fc + 1 and the elapsed time is t-(fc) + 1. 
Thus T^^^' obeys the recursion 

= E+ik)U{k) + E^ik) [Uk) + 1 + ri^^)] . (40) 

If we allowed the walk to hop to the left (with probability 
q) when the leftmost cookie is reached, then Eq. H4()|l 
would be replaced by the exact equation 

= i?+(fc)t+(fc) + i?_(fc)[t_(fc)+p(l + Tig) 

+ + + 9^3 + Ti°^1) + ...]• 

The additional terms represent the contribution of paths 
in which the random walk makes 1, 2, 3, etc. consecutive 
steps to the left upon hitting the last cookie before a step 
to the right is made. This exact equation is not recursive 
and for simplicity we therefore consider the approxima- 
tion recursion for the lower bound in what follows. 
By iterating (|^ . we express T^^^^ compactly as 

oo e 

T^'^^'^-r, + ^rfe+,ni?_(fc + j-l) 

£=1 j=l 

oo e 

£=1 j = l 

where Tk = E^{k)t^{k) + E+{k)t+{k) ^ 2(fc - 2) is the 
unconditional exit time, namely, the time for the ran- 
dom walk to exit either side of the interval [0, fc] when 
the starting point is x = 2. Using this result for r^, 
as well as that for E-{k) in Eq. (|41() . we obtain after 
straightforward algebra: 

Tj:^^^ = 3(fc - 2) + 2(fc - l)(fc - 2) ^ - . (42) 

j=k ^ 



The upper limit in the sum accounts for case where a 
cookie-free defect of length fc occurs close to the right 
edge of the domain [0, L] so that terms or order L are 
needed to account for the growth of the cookie-free de- 
fect. Thus we infer that mean time for a 2-excited walk 
to traverse a finite system of length L scales at least as 
fast as LlnL. This result is then consistent with {x{t)) 
growing slower than linearly with t. 



VI. SUMMARY 

The excited random walk is a very rich stochastic pro- 
cess in which the random walk motion is influenced by 
a transitory bias. Each site of a lattice is initially pop- 
ulated with a finite number of cookies. When a walk 
reaches a site that contains cookies, the walk eats the 
cookie and there is a p/q right/left bias at the next step. 
On the other hand, the walk undergoes unbiased motion 
if the visited site does not contain cookies. Thus as the 
cookies are eaten, the bias gradually disappears. If a 
walk is recurrent, then the interplay between the recur- 
rence and the bias is particular delicate and the excited 
random walk will have very different behavior from pure 
diffusion or from diffusion with a constant bias. For this 
reason, the excited random walk appears to be most in- 
teresting in one dimension. 

In the simplest case of the 1-excited walk on the posi- 
tive half-line, we employed a physical approach, comple- 
mentary to that of Dickman et al. [lOl [ll| , to show that 
the first-passage probability for the walk to hit the ori- 
gin at time t decays as t~^~'^ for any q > 0. Thus the 
1-excited walk is always recurrent, except for the trivial 
case of perfect bias p — 1- There is especially paradoxical 
behavior in the case ofp < 1/2, where the cookie-induced 
bias is directed toward the origin. Here the walk exhibits 
two seemingly incompatible features: the mean time for 
the walk to hit the origin is infinite, while the walk strays 
only a finite distance from the origin during this infinite 
excursion. 

On the infinite line, the 1-excited walk has the same 
first-passage characteristics as on the semi-infinite line. 
However, the probability distribution of the walk and 
the probability distributions for the positions of the last 
cookie exhibit unusual features. For the probability dis- 
tribution of the walk, there appears to be a singularity at 
the origin that does not obey scaling while the rest of the 
distribution does scale. There is also an asymmetry in 
the span of the walk that we were able to characterize by 
a probabilistic argument. The probability distributions 
for the positions of the leftmost and rightmost cookies 
have power law tails that we have been unable, thus far, 
to explain. These tails seem to be connected with the 
power-law behavior for the first-passage probability it- 
self. 

The 2-excited walk shows a wide range of strange phe- 
nomena. A transition between (conventional) recurrence 
and transience occurs as p passes through a threshold 
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value Pc = 3/4. For p < pc, the 2-excited walk ap- 
pears to qualitatively resemble the 1-excited walk. On 
the other hand, for p > pc, the 2-excited walk is tran- 
sient and this attribute appears responsible for several 
unexpected features, such as the lack of scaling in the 
probability distribution of the walk and in the spatial 
distribution of cookies, as well as the existence of very 
different distributions of 1-cookie and 0-cookie domains. 
Most of these features do not yet have clean explanations 
and the model is ripe for further exploration. 



sufficiently close to the position of the last cookie: 
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APPENDIX: MASTER EQUATION 

In this appendix, we outline a Master equation ap- 
proach to determine the survival probability of the 1- 
excited random walk on the infinite half-line. We me rely 
follow the method introduced by Dickman et al. [lO, 
to describe an unbiased walk between a fixed absorbing 
boundary at a; = and a moving reflecting boundary 
that moves to the right with a specifled probability each 
time it is hit by the walk. This perturbed random walk 
motion is essentially the same mechanism as that of the 
excited random walk. While the excited walk is ostensi- 
bly not Markovian, the process can be made Markovian 
by considering an extended state space that is defined 
by the position of the walk x and also the positions of 
the leftmost and rightmost cookies, y and y. The ex- 
act master equation for the joint probability distribution 
of the position of the walk and the positions of the last 
cookies P{x,y, y, t) can be obtained in terms of a Markov 
process. 

We consider the simpler problem of the 1-excited ran- 
dom walk with an absorbing boundary at the origin. The 
system is now described by the joint probability for the 
cookie to be at x and the leftmost cookie to be at y, 
P{x,y,t). This distribution obeys the recursion formula 

Pix, y,t+l) = ^P{x - 1, y, t) + ^P{x + 1, y, t) , (A.l) 

for a; = l,2,...,?/— 3. At the origin there is an absorbing 
boundary condition 



P(0,y,0 = , 



(A.2) 



while the effect of the cookies is described by slightly 
modified equations for the probability when the walk is 



These equations are defined over the wedge domain x > 
0,x < y. We choose the initial conditions that the walk 
starts a.t X = 1 and that the last cookie is also at a; = 
1; thus P{x,y,0) — S^^iSy^i. These equations can be 
iterated then numerically and this is the method that we 
used for the simulations presented in Sec. II. 

From these Master equations, we compute the proba- 
bility for the excited random walk to reach the origin for 
the first time at time t, from which the survival proba- 
bility of the walk follows. We merely outline this deriva- 
tion as it very similar to that given by Dickman et al. 
[Tol ITU . We first rewrite the Master equations in terms 
of the generating function P(x, y, z) = X^t^o z*P{^i 2/i 0- 
In this representation, the solution can be written in the 
product form 



P{x,y,z) = A{x,z)B{y,z), 



(A.4) 



where 



i(a;, z) = - A"^ , with A = ^ + \/ 4 - 1, (A.5) 

Z V 2 



and where B is defined recursively by 



B{y + l,z) _ tanh Ay -I- 2gf A 



B{y.z) 



tanh Ay + A 



(A.6) 



in the e — z — \ ^ Q limit, with A = In A. 

The long-time behavior of the survival probability 
^{^) ^ X^a; y P{^^ 2/) ^) ^au uow bc obtalncd from the 
z — ^ 1 limit (A 0) of P{x,y, z). For any finite y, and 
thus for y = yo ^ 2(p - 1), B{yo) B{1) ^ 1/A as 
A ^ 0. On the other hand, for y > yo we approximate 
the recursion formula ljA.6|l by 



In 



Bjy) 

Biyo) 



El- 

k=Vo 



l + 2{p- l)A/tanhAfc 
1 + A/tanhAA: 



(2p-3)A^ 



1 



tanh Ak 



(A.7) 



(3 - 2p) In 



k=Vo 

sinh Ayo 
sinh Ay 
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Putting these terms together we obtain 

oo y 
Siz) = ^i3(y,z)^i(a;,z) 

y—1 x—0 

B(y,z)isinh2^ 



E 



A1-2P J2 



sinh^ Ay/2 
^ sinh^-^f Ay 



A-2p 



sinh^ u/2 



Ayo 



sinh 



3-2p, 



Since the integral above goes to a finite constant as A 
0, S{z) ~ A-2p - e-P, which corresponds to S{t) ~ 



For the 1-excited random walk on the infinite line the 
(A. 8) phase space is now the three-dimensional domain defined 
by y < x < y. In principle, we can write recursion for- 
mulae similar to Eqs. (|A.ip - (|A.3|I for P{x,y,y,t), but 
the solution to these equations no longer seems practical. 
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